Measurement of the Rate of Muon Capture in Hydrogen Gas and 
Determination of the Proton's Pseudoscalar Coupling g p 
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The rate of nuclear muon capture by the proton has been measured using a new experimental 
technique based on a time projection chamber operating in ultra-clean, deuterium-depleted hydrogen 
gas at 1 MPa pressure. The capture rate was obtained from the difference between the measured 
yU~ disappearance rate in hydrogen and the world average for the /i + decay rate. The target's low 
gas density of 1% compared to liquid hydrogen is key to avoiding uncertainties that arise from the 
formation of muonic molecules. The capture rate from the hyperfine singlet ground state of the p,p 
atom is measured to be A s = 725.0 ± 17.4 s _1 , from which the induced pseudoscalar coupling of 
the nucleon, g P (q 2 = — 0.88 m 2 ,) = 7.3 ± 1.1, is extracted. This result is consistent with theoretical 
predictions for g p that are based on the approximate chiral symmetry of QCD. 

PACS numbers: 23.40.-s, 11.40. Ha, 13.60.-r, 14.20. Dh, 24.80. +y, 29.40.Gx 



We report the first result of the MuCap experiment for 
the rate A s of the semileptonic weak process of ordinary 
muon capture (OMC) by the proton, 

H~~ + p — ► n + . (1) 

This fundamental process, like neutron beta decay, in- 
volves the vector and axial- vector form factors g v {q 2 ) 
and g A {q 2 ), which characterize the microscopic QCD 
structure of the nucleon in electroweak charged-current 
interactions. Due to its larger momentum transfer q 2 = 
—0.88 m 2 , reaction (fil is also sensitive to the weak mag- 
netic and pseudoscalar induced form factors, g M (q 2 ) and 
g P {q 2 ). Form factors g y (<? 2 ), g M (Vo) and Sa{<Iq) are 
accurately determined by experimental data and Stan- 
dard Model symmetries and contribute an uncertainty of 
only 0.46% to A s pQ. Process ([!]) provides the most di- 
rect probe of g P = gp(g 2 ), the pseudoscalar coupling of 
the nucleon's axial current, which is by far the least well 
known of these form factors. 

The form factor g P {q 2 ) arises mainly from the coupling 
of the weak leptonic current to the nucleon via an inter- 
mediate pion, which generates a pole term that domi- 
nates at (?g. Early theoretical expressions for g P were 
derived using current algebra techniques; now g P can 
be systematically calculated within heavy baryon chiral 
perturbation theory (HBChPT) up to two-loop order [2]. 
The precise QCD result g p = 8.26 ± 0.23 [1 i] follows 




^„ P (io 3 s- 1 ) 



FIG. 1: Experimental and theoretical determinations of g p , 
presented vs. the ortho-para transition rate A op in the p\ip 
molecule. The most precise previous OMC experiment [7] 
and the RMC experiment [8] both depend significantly on 
the value of A op , which itself is poorly known due to mutually 
inconsistent experimental (A^* 1 [5], A^p 2 [TO]) and theoretical 
(Xjp El) results. In contrast, the MuCap result for g P is 
nearly independent of molecular effects. 

from the basic concepts of explicit and spontaneous chiral 
symmetry breaking, and thus its experimental confirma- 
tion is an important test of QCD symmetries jlHB]. 

Experimental OMC efforts span a period of more 
than forty years, and more recently radiative muon cap- 
ture (RMC) by the proton was measured for the first 
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time [Sj. However, as shown in Fig. [I] the situation prior 
to the present experiment was inconclusive, as the re- 
sults lacked sufficient precision due to ambiguities in the 
interpretation as well as technical challenges. 

The problems of interpretation can be appreciated by 
considering the chain of reactions possible for negative 
muons after stopping in a hydrogen target of density <f> 
relative to liquid hydrogen (LH 2 ) [5]. Stopped muons 
immediately form ground state fip atoms whose hyper- 
fine states are populated in a statistical manner. The 
upper triplet spin state is rapidly depopulated in col- 
lisions with H 2 molecules, and for densities <f> > 0.01 
all muons reach the ixp singlet state well before 100 ns. 
From there, muons can either decay with a rate close 
to A+ = J_/t+ w 0.455 x 10 6 s _1 , or be captured via 
reaction (fil at the predicted rate A s « 710 s _1 . Com- 
plications arise at higher densities, however, as \xp atoms 
increasingly collide with target H 2 molecules to form p/j,p 
molecules. The p\xp molecules arc initially created in the 
ortho state at the density-dependent rate </>A f, where 
A f 2.3 x 10 6 s _1 , and then de-excite to the para state 
at rate A op . The nuclear capture rates from the ortho and 
para states, A om » 506 s _1 and A pm ss 208 s _1 [4], are 
quite different from each other and from A s , so knowl- 
edge of the relative populations of the \xp and pfip states 
under any particular set of experimental conditions is 
crucial for a correct determination of g P . Alas, A op is 
poorly known [9TI11]. This prevents a clear interpreta- 
tion of the most precise OMC experiment [7 a , which was 
performed in LH 2 where muon capture occurs predomi- 
nantly in pfip molecules. The RMC process is less sensi- 
tive to A op , but the large molecular uncertainties make it 
difficult to draw firm conclusions from the RMC experi- 
ment [8], whose results initially suggested a nearly 50% 
higher value for g p than predicted. 

Direct measurement of A s is technically difficult be- 
cause process ([T|) is rare (branching ratio = 0.16%) and 
leads to an all-neutral final state. Moreover, target impu- 
rities and muon stops in detector walls must be scrupu- 
lously avoided, as negative muons preferentially and ir- 
reversibly transfer from /ip to heavier elements, and the 
nuclear muon capture rate increases roughly proportional 
to Z 4 . The two previous muon capture experiments using 
low-density gas targets and neutron detectors obtained a 
precision in A s of 9% [T2] and 13% [13], respectively. 

The MuCap experiment employs novel techniques to 
minimize or avoid many of the problems described above. 
The measurement is performed using hydrogen at den- 
sity <fi — (1-12 ± 0.01) x 10~ 2 , where pfip formation is 
slow and 96% of all captures proceed from the \±p singlet 
state. The significant background from muon stops in 
wall materials, inherent when using a low-density target, 
is eliminated by reconstructing the muon stopping point 
in an active target consisting of a hydrogen time projec- 
tion chamber (TPC). The capture rate is determined us- 
ing the lifetime technique [7], that is, from the difference 




FIG. 2: Simplified cross-sectional diagram of the MuCap de- 
tector. The detector components are described in the text. 

between the measured disappearance rate A~ « A+ + A s 
of negative muons in hydrogen and the fi + decay rate A+ , 
where it is assumed that free /i~ and ii + decay with iden- 
tical rates according to the CPT theorem. 

The experiment was conducted in the 7rE3 beamlinc 
at the Paul Scherrer Institute, using a ss 20 kHz DC 
muon beam tuned to a central momentum of 32.6 MeV/c. 
As illustrated in Fig. [2] incident muons first traverse a 
plastic scintillator (/iSC) and a multiwire proportional 
chamber (/iPC), and then pass through a 0.5-mm-thick 
hemispherical beryllium window to enter an aluminum 
pressure vessel filled with ultra-pure, deuterium-depleted 
hydrogen gas at a pressure of 1.00 MPa and at ambient 
room temperature. In the center of the vessel is the TPC 
(sensitive volume 15 x 12 x 28 cm 3 ), which tracks incom- 
ing muon trajectories and thus enables the selection of 
muons that stop in the gas at least 15 mm away from 
chamber materials. Approximately 65% of the muons 
passing through the [iSG stop within this fiducial volume. 
The ionization electrons produced by incoming muons 
drift downwards at velocity 5.5 mm//is in an applied field 
of 2 kV/cm, towards a multiwire proportional chamber 
containing perpendicular anode and cathode wires. The 
anode plane consists of wires with 25 fj.ni diameter and 
4 mm spacing, and a high voltage of 5.0 kV across the 
3.5 mm half-gaps achieves a moderate gain of 60 in hydro- 
gen. Digital signals from three-level discriminators are 
recorded, with the energy thresholds adjusted to trigger 
on (i) fast muons, (ii) the Bragg peaks near the muon 
stopping points, and (iii) the larger energies that may be 
deposited by recoiling nuclei following muon capture by 
gas impurities. 

The TPC is surrounded by two cylindrical wire cham- 
bers (ePCl, ePC2), each containing anodes and in- 
ner/outer cathode strips, and by a hodoscope bar- 
rel (eSC) consisting of 16 segments with two layers of 
5-mm-thick plastic scintillator. This tracking system de- 
tects outgoing decay electrons with 37r solid angle ac- 
ceptance. All data are recorded in a triggerless, quasi- 
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a. No u-e impact cut 

b. 120 mm u-e impact cut 

c. 50 mm u-e impact cut 
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FIG. 3: Lifetime spectra of negative muons. The signal-to- 
background ratio improves with tighter cuts on the p,-e vertex. 

continuous mode to avoid deadtime distortions to the 
lifetime spectra. Custom-built time-to-digital convert- 
ers (TDCs) digitize hit times for the TPC and the elec- 
tron wire chambers. The muon and electron times 
and t e are established by the /iSC and eSC detectors, 
and recorded in separate CAEN V767 TDC modules. 

All TPC materials were carefully selected for high vac- 
uum operation. Prior to the run, the TPC system was 
heated to 115°C under vacuum for several weeks to re- 
move impurities. The system was filled with deuterium- 
depleted hydrogen through a palladium filter to remove 
impurities. During data taking, the gas was continu- 
ously circulated via an adsorption cryopump system and 
cleaned by cooled Zeolite filters [T3], which achieved an 
equilibrium concentration (by number) of c z < 5 x 10~ 8 , 
as monitored by direct TPC detection of recoil nuclei 
from muon capture by impurities. Gas chromatography 
measurements established that the atomic concentration 
of nitrogen was below 10 -8 , and the post-run installa- 
tion of a humidity sensor with 10 -9 sensitivity into the 
gas circuit indicated that the primary contaminant was 
H 2 outgassing from within the pressure vessel. 

The isotopic purity of the hydrogen is critical. Muons 
preferentially transfer from /ip to /id at the rate <j)Cd\ P d, 
where ca is the deuterium concentration and X p( i w 
1.4 x 10 10 s _1 . Whereas fip diffusion is on the order 
of mm, /id atoms can diffuse cm-scale distances due to 
a Ramsauer-Townsend minimum in the fid + p elastic 
scattering cross-section. As a result, fj,d atoms can drift 
sufficiently far away from the muon's original stopping 
point that the decay event will be rejected by the fi-e 
vertex reconstruction cut in a time-dependent manner. 
In addition, \id atoms can drift into surrounding materi- 
als and be captured there. Our target gas was produced 
via electrolysis of deuterium-depleted water, and accel- 
erator mass spectrometry (AMS) measurements |15j de- 
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Z > 1 impurities 
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5.0 


[t,d diffusion 


-10.2 


1.6 


[ip diffusion 


-2.7 
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Total 


-32.1 


8.5 



TABLE I: Systematic corrections and uncertainties applied to 
the observed /x~ disappearance rate A. 

termined that c d = (1.44 ± 0.13) x 1CT 6 , roughly 100 
times below deuterium's natural abundance. This result 
was independently confirmed from our data, by analyzing 
the observed losses of muon decay events as a function of 
the imposed fi-e vertex cut. 

The time differences between muon arrivals and de- 
cay electron emissions, At = t e — t^, are histogrammed 
into lifetime spectra (Fig. |3j. Only muons that are sep- 
arated in time by ±25 ps from other muon arrivals are 
accepted. While this condition cuts the usable statistics 
by ~ 68%, it is essential for avoiding systematic distor- 
tions to the background which can arise from ambigui- 
ties in resolving multiple muon tracks in the TPC, and 
it dramatically improves the signal-to-background ratio. 
As shown in Fig. [3] further background suppression can 
be achieved by performing a vertex cut on the impact 
parameter between each decay electron's trajectory and 
its parent muon's stopping point. In the final analysis 
we employ a loose impact parameter cut of 120 mm as 
an optimal compromise between the competing demands 
for a good signal-to-background ratio and minimization 
of losses due to [id diffusion out of the cut volume. 

We fit the fi~ lifetime spectra with the simple expo- 
nential function f(t) = Nw\e~ xt + B, where the free pa- 
rameters are the number of reconstructed decay events 
N, the disappearance rate A, and the accidental back- 
ground level B\ w is the fixed 40 ns histogram bin width. 
We studied an assortment of analysis conditions, includ- 
ing different time ranges (0.1-24 /is is typical), fiducial 
cuts, and detector combinations, and typically obtained 
X 2 /dof = 0.95-1.02 for 600 degrees of freedom. 

In reality, the experimental pT lifetime spectrum is not 
a pure exponential, but has a more complicated shape 
due to contributions from p[ip molecules and hydrogen 
gas impurities. However, these effects are sufficiently 
small that their perturbations A A to the exponential de- 
cay rate A are linear and can be corrected sequentially. 
The main corrections to A were derived directly from ex- 
perimental data, with some additional information from 
external measurements and literature. For residual c z 
below a few times 10 -6 , the correction AA Z scales with 
the observed impurity capture yield per muon, Y, 
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yield Y z w 11 x 10 -6 , received contributions from nitro- 
gen and humidity in weights of approximately a N = 0.05 
and a = 0.95. The factors (^A) callb were empirically 
fixed by calibration runs involving N 2 and O (in the form 
of H 2 0) concentrations 50-1000 times above their values 
in the clean fill. We find AX Z = —19.2 ± 5.0 s -1 , where 
the error is dominated by a conservative estimate of the 
^AA-jcahb vame f or q determined during our 2006 run- 
ning period. The correction for deuterium-related diffu- 
sion effects, AXd = —10.2 ± 1.6 s _1 , was obtained by a 
zero-extrapolation procedure using data from a run with 
a hydrogen filling of c d = (122 ± 5) x 10~ 6 . 

The preceding corrections are summarized in Table [i] 
There we also present four additional sources of uncer- 
tainty, including a conservative error of 5 s _1 that ac- 
counts for the spread in results observed for a variety 
of consistency studies, as performed by two independent 
analyses. To prevent bias, the master clock (accurate to 
10 -8 [16]), was detuned by an offset, which was concealed 
until the data analysis was complete. 

The final result for the yT disappearance rate in pure 
hydrogen, based on N = 1.6 x 10 9 fully tracked, pileup- 
protected decay events from our 2004 data set, is A~ = 
455 851.4 ± 12.5 sta t ± 8.5 syst s _1 . As a consistency check, 
we also measured the \i + decay rate from N = 0.5 x 10 9 
events to be A+ = 455 164 ± 28 s _1 , in agreement with 
the world average. 

The observed pT disappearance rate can be written as 

K = (A+ + AA W ) + A s + AA PW . (2) 

Here AA W = —12.3 s _1 describes a small reduction in 
the muon decay rate in the bound /ip system [T? . The 
term AA PM , = — 23.5±4.3±3.9 s _1 accounts for captures 
from pfip molecules, and is calculated from the full 
kinetics in pure hydrogen. Its error terms come from 
our estimates A f = (2.3 ± 0.5) x 10 6 s _1 and A op = 
(6.9±4.3) x 10 4 s _1 , respectively, which cover most of the 
existing literature values. As muon capture from the fj,p 
singlet component dominates both in pp atoms and ppp 
molecules, AA ppp implicitly depends on As, which leads 
to a 3.2% loss in sensitivity when determining As from 
Eq. pj. Using the new world average A+ = 455 162.2 ± 
4.4 s [16j . we determine the rate of muon capture by 
the proton to be 

A MuCap = ?25 Q ± ^ ± 10 Jsyst g -l (3) 

To compare with theory we consider the two recent 
NNLO calculations of A s , 687.4 s" 1 [IB] and 695 s" 1 [10], 
here averaged to 691.2 s _1 . Adding the very recently 
calculated radiative correction A# = 19.4 s _1 [20J (in- 
creased from Aji = 4.5 s" 1 [2T]) yields the value A^ 11 = 
710.6 s _1 and enables us to calculate 

g MuCap = gTh + ||P ^MuCap _ ^ = 7 .3±1.1 , (4) 



where gj^ = 8.26 0], ^ = -0.065 s [5J, and only 
the experimental uncertainty from Eq. ([3]) is propagated. 
The linear expansion in Eq. Q is valid because of the 
small difference Ag U ap — A™, but should be refined 
once further theoretical work clarifies the present 1% dif- 
ference between calculations [T5J and [TO] and quantifies 
all sources of theoretical uncertainty at the sub-percent 
level. 

The current information on g p is summarized in Fig. [TJ 
the constraints [5] from the OMC experiment [7] are up- 
dated to reflect the larger An [20]. The situation before 
MuCap was inconclusive and exhibited mutually incon- 
sistent theoretical predictions and experimental determi- 
nations of both g p and A op . The low gas density in Mu- 
Cap renders our result relatively insensitive to A op and 
thus avoids most model dependence, enabling us to re- 
port the first precise, unambiguous determination of g P . 
This experimental result agrees with present theory to 
within lcr and does not support a dramatic deviation 
from the chiral prediction as the RMC result originally 
had implied. Additional data are being collected with 
the aim of a more than twofold reduction of statistical 
and systematic uncertainties. 
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